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摘要 : 谱 线 观测 在 天 文 领域 起 着 至 关 重 要 的 作用 ， 较 高 的 频谱 分 辩 浴 可 以 为 研究 天 体 细 
节 带 来 便利 ,但 在 高 频段 谱 线 观测 领域 ,使 用 传统 的 采样 方法 实现 较 高 的 谱 线 分 辨 率 十 分 
难 。 目 前 ， 一 种 利用 高 性 能 芯片 实现 的 用 于 射电 天 文 厘米 -分 米 波 谱 线 观测 领域 的 新 方法 克 
服 了 传统 方法 的 诸多 弊端 ， 使 用 较 低 的 采样 带宽 就 可 实现 较 高 的 频谱 分 辨 率 。 这 是 一 套 完整 
的 采集 、 处 理 、 存 储 方案 。 首 先 将 高 频 信 号 进行 模拟 下 变频 处 理 ， 再 将 处 理 过 的 低频 信号 交 
由 上 位 机 软件 进行 快速 全 里 叶 变换 、 积 分 处 理 ， 最 后 保存 成 需要 的 文件 格式 。 整 个 系统 的 采 
集 参 数 由 观测 者 通过 上 位 机 软件 进行 配置 。 
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谱 线 观测 在 获得 天 体 信息 方面 具有 重要 意义 。 对 于 高 频段 的 观测 要 求 数据 采集 设备 具有 更 高 的 采 
样 速率 。 各 大 芯片 公司 也 在 研制 更 高 速 更 高 带宽 的 数据 采集 芯片 ， 其 中 一 部 分 世 斤 已 经 应 用 于 射电 天 
文 观测 。 然 而 ， 高 速 数据 采集 带 来 的 头 端 使 得 资源 浪费 在 一 些 没有 可 观测 信号 的 频率 上 ， 而 且 采 样 率 
越 高 ， 浪 费 越 多 。 在 集成 电路 日 新 月 异 的 今天 ， 一 种 新 的 观测 控制 方法 应 运 而 生 IT 


1 技术 背景 


谱 线 在 射电 天 文 观测 中 具有 举足轻重 的 地 位 ， 可 用 于 诊断 相关 天 体 的 基本 物理 及 化 学 组 成 ， 如 激 
发 温度 、 粒 子 密 度 、 物 质 运动 区 的 速度 场 、 磁 感应 强度 及 各 种 元 素 的 化 学 丰 度 等 。 通 过 观测 获得 分 子 
云 及 恒星 形成 区 的 谱 线 数 据 后 ， 可 以 通过 适当 的 近似 方法 求解 统计 平衡 辐射 转移 方程 ， 从 而 获得 天 体 
的 各 种 信息 。 

通过 对 分 子 谱 线 的 观测 ， 能 够 研究 包括 恒星 形成 区 早期 、 恒 星 演化 晚期 、 超 新 星 爆炸 (超新星 遗 
迹 ) 、 行 星 状 星云 、 芷 星 、 河 外 星系 等 一 系列 的 天 文学 现象 。 对 它们 在 C 波段 谱 线 的 观测 研究 ， 关 联 
了 很 多 研究 内 容 ”1 ， 例 如 : (1) 基 于 6.7 GHz 甲醇 脉 泽 相关 的 研究 ，6.7 GHz 甲醇 脉 泽 由 于 其 辐射 流 
量 大 、 致 密 ( 角 大 小 约 4 毫 角 秒 ， 线 大 小 约 几 个 天 文 单位 ) 、 寿 命 长 ( 约 4 年 )、 普 遍 存 在 (银河 系 内 探 
测 到 1 000 多 个 源 ) 以 及 非常 小 的 内 部 固有 运动 等 特点 成 为 天 文学 和 天 体 测 量 学 中 的 一 大 热点 。 它 们 
不 仅 被 用 来 开展 与 大 质量 恒星 形成 有 关 的 一 系列 人 研究， 而 且 通 过 对 它们 视差 的 测量 可 以 精确 确定 银河 
系 的 一 些 关 键 参 数 ( 如 旋 臂 结构 ， 这 是 迄今 为 止 最 高 精度 的 测量 方法 ) ， 是 银河 系 研究 的 重点 ; (2) 基 
于 产 基 和 甲醛 谱 线 的 研究 也 是 近来 受到 广泛 关注 的 热点 领域 ， 包 含 甲 醛 样本 搜寻 以 及 他 们 的 辐射 特 
性 、 产 基 脉 泽 的 搜寻 以 及 通过 产 基 偏振 观测 研究 磁场 等 内 容 ; (3) 另 外 还 有 很 多 研究 内 容 ， 比 如 C X 
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段 的 分 子 谱 线 搜寻 等 。 

高 频谱 分 辨 率 的 谱 线 观测 可 用 于 研究 目标 天 体内 部 更 细节 的 信息 ， 它 们 通常 可 以 追踪 到 更 小 尺度 

4 结构 以 及 更 详细 的 动力 学 特征 等 。 在 某 些 高 精度 的 谱 线 观测 中 需要 极 高 的 频率 分 辨 率 ， 例 如 对 塞 曼 
Ps 的 观测 中 ， 采 用 中 性 氢 HI 观测 ， 每 微 高 斯 磁场 产生 的 频率 移动 2.8 Hz; 如 采用 OH 的 基态 ， 每 
微 高 斯 磁场 在 1.665 GHz 谱 线 上 产生 的 频率 移动 只 有 3.27 Hz, MÆ 1.667 GHz 谱 线 上 产生 的 频率 移 
动 是 1. 96 Hz。 然 而 ， 由 于 传统 的 射电 天 文 终端 设备 是 对 射电 天 文 模拟 信号 的 高 速 模 拟 数字 转换 需 过 
采样 原理 实现 : 采样 率 为 人 的 模拟 数字 转换 希 大 约 能 上 覆盖 人 .2 带宽 的 信号 ， 目 前 增加 采样 率 成 
为 解决 大 市 宽 数 据 采 集 、 en NI E UN, 

根据 多 普 勒 关系 ， 对 于 给 定 的 速度 弥散 Av， 对 应 的 频谱 宽度 应 为 


EET (1) 
根据 采样 率 信 ， 、 快 速 傅 里 叶 变换 点 数 N 和 频谱 分 辨 率 Af 三 者 的 关系 中， 
As (2) 


速度 分 辨 率 与 采样 率 、 人 对 同一 条 谱 线 在 同样 入 点 
数 情况 下 ， 采 样 率 越 高 ， 速 度 分 辨 率 越 低 “” : 


人 le 
Av =, (3) 
Nf, 


在 传统 采样 原理 下 ， 谱 线 的 高 精度 轮廓 结构 观测 非常 困难 ， 因 为 在 采样 率 一 定 的 情况 下 ， 要 得 到 
更 高 的 频谱 分 辨 率 就 必须 增加 快速 傅 里 时 变 换 点 数 W， 这 就 需要 耗费 处 理 需 内 部 乘法 禹 、 累 加 器 以 及 
存储 单元 等 资源 ， 大 大 增加 了 宛 余 数据 量 ， 对 无 信号 频段 也 需要 同时 进行 分 析 。 

目前 ， 各 大 芯片 公司 和 模块 化 仪器 公司 相继 进行 更 高 速 更 高 带宽 的 数据 采集 芯片 和 设备 设计 开发 ， 
多 种 高 速 数 据 采集 设备 已 经 应 用 于 天 文 观测 。 例 如 : Acqiris 公司 开发 的 高 速 数据 采集 卡 AC240, Honivu 
为 1 GHz@8 bit， 已 经 应 用 于 多 台 射 电 望 远 镜 ， 包 括 德 令 哈 的 13.7 m ZKE, KOSMA ZKK V. 


毫米 波 望 远 镜 ; B. Klien 等 人 研发 的 FFTS :5 数据 采集 终端 均 是 目前 采样 速度 最 快 的 射电 天 文 终端 。 
根据 原子 受 激 辐射 的 特点 ， 分 子 谱 线 呈 现 离 

散 局 部 密集 线 谱 的 特点 ， 在 其 他 频段 能 够 观测 的 NEN 

射电 天 文 谱 线 又 非常 少 。 这 样 ， 高 速 数据 采集 的 autocorrelation a 

优势 就 浪费 在 一 些 没有 可 观测 信和 号 的 频率 上 ， 且 25 EE | 


XFFTS: 2.5 GHz/64K 
KERRE, WARS, 

2008 年 文 [4] 实 现 的 100 MHz 输入 带宽 的 数 
字 终 端 TasPGA ^ ， 观 测 脉 泽 源 G285. 35 +0. 00 
在 8 MHz 带宽 内 实现 488. 3 Hz 的 频谱 分 辨 率 需 
要 16 384 点 快速 傅 里 叶 变 换 ; 文 [5] 在 FFTS 上 
实现 88.3 KHz， 需 要 32 768 点 快速 傅 里 叶 变 换 ; 
文 [5] 总 结 的 随 着 采集 速度 的 上 升 ， 达 到 同样 
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针对 上 述 问题 ，2005 年 文 [6] 尝试 在 埃 费 D = | 
尔 斯 贝 格 ( Effelsberg) 100 m 射电 望远镜 上 建立 2. 4 8 16 64 128 
基于 DDC-GCA016 数 学 TRE is 频 p ( Digital Dow Number of spectral channels [x 1024] 
Converters, DDC ) 的 数字 频谱 分 析 终 端 ICS 554C, 图 1 快速 傅 里 叶 变 换 点 数 关系 图 
引入 了 直接 数字 下 变频 的 概念 与 谱 线 数字 终 Fig.1 The relationship between bandwidth 


5] 如 图 2. 它 有 两 种 工作 方式 ”9 ; (1) 全 带 and FFT points 
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WRIN: 50 MHz， 直 接 由 模拟 数字 转换 器 采集 进入 现场 可 编程 门 阵列 做 快速 依 里 叶 变 换 分 析 ; (2)78 
带 工作 方式 ， 经 数字 下 变频 器 选 通 20 KHz-10 MHz 带宽 ， 再 送 入 现场 可 编程 门 阵列 进行 处 理 。 两 种 
方式 对 同一 条 谱 线 的 观测 效果 如 图 3、 图 4。 
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图 2 数字 频谱 分 析 终 端 ICS 554C 
Fig.2 Digital spectrum analysis terminal ICS 554C 
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图 3 全 带宽 模式 下 对 中 性 所 和 卷 基 脉 泽 源 的 观测 
Fig.3 Observation of Neutral hydrogen and OH in full bandwidth mode 
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图 4 RPA TFE 2 MHz 带宽 对 中 性 氨 和 券 基 脉 泽 源 的 观测 
Fig.4 Observation of Neutral hydrogen and OH with a 2MHz bandwidth in narrow band mode 
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在 2 MHz WERA F, HATEEN BR PII es fS E p AERÉd X dE I 25 倍 的 频率 分 
辨 率 。 同 时 ， 文 [7-8] 提 出 了 基于 欠 采 样 方案 的 采集 终端 ， 利 用 模拟 数字 转换 器 的 欠 采 样 性 可 以 达到 
更 高 的 分 辨 率 。 但 是 以 上 两 种 方案 都 以 增加 电路 的 复杂 结构 为 代价 ， 在 集成 电路 日 新 月 异 的 今天 ， 新 
的 集成 电路 能 够 取代 上 述 方案 2。 

目前 ADI 公司 推出 的 AD93xx 系列 芯片 是 一 种 高 性 能 、 高 度 集成 的 RF Agile Transceiver" 捷 变 收 
发 需 。 器 件 的 可 编程 性 和 宽带 使 其 成 为 多 种 收发 器 应 用 的 理想 选择 。 需 件 集 射频 前 端 与 灵活 的 混合 信 
号 基带 为 一 体 ， 集 成 频率 合成 器 ， 为 处 理 需 提 供 可 配置 数字 接口 ， 从 而 简化 设计 。 例 如 : AD9361 T. 
作 频 率 为 70 MHz 至 6.0 GHz， 涵 盖 大 部 分 特许 执照 和 免 执 照 频 段 ， 支 持 的 通道 带宽 从 不 到 200 KHz 至 
56 MHz", 

较 之 于 传统 的 高 速 采 样 数字 频谱 终端 ， 根 据 奈 奎 斯 特定 律 达到 6 GHz 的 输入 带宽 至 少 要 有 12 Gsps 
的 采样 速度 ， 为 此 在 捷 变 收发 器 选择 200 kHz 输入 带宽 的 情况 下 ， 到 达 同 样 的 分 辨 紊 ， 相 当 于 点 数 提 
高 了 240 000 售 。 

系统 采用 先进 的 集成 射频 捷 变 收发 右 人 研发 新 型 的 专用 谱 线 接收 机 ， 实 现 快速 宽带 内 可 切换 的 射电 
天 文 终端 及 相关 软件 ， 并 初步 为 云南 天 文 台 40 m 射电 望远镜 提供 高 精度 谱 线 观测 。 


2 工作 原理 | 


RXIA P. 


2.1 捷 变 收发 硕 AD9361 

AD9361 是 捷 变 收发 方案 的 典型 代表 ， 其 工 
作 频 率 为 70 MHz 至 6.0 GHz, 支持 的 通道 带宽 
从 不 到 200 kHz 至 56 MHz。 信 号 处 理 板 如 图 5。 

两 个 独立 的 直接 变频 接收 器 拥有 首屈一指 的 
噪声 系数 和 线性 度 。 每 个 接收 子 系统 都 有 独立 的 
自动 增益 控制 (Automatic Gain. Control, AGC). 
直流 失调 校正 、 正 交 校 正和 数字 滤波 功能 ， 从 而 
消除 了 在 数字 基带 中 提供 这 些 功能 的 必要 性 。 
AD9361 还 拥有 灵活 的 手动 增益 模式 ， 文 持 外 部 
控制 。 每 个 通道 搭载 两 个 高 动态 范围 的 模拟 数字 
转换 器 ， 先 将 收 到 的 工 信 号 和 Q 信号 进行 数字 化 
处 理 ， 然 后 将 其 传 给 可 配置 抽取 滤波 器 和 128 抽 TXB NY” 


DATA INTERFACE 


RADIO 
(SWITCHING 


头 有 限 脉 冲 啊 应 (Finite Impulse Response, FIR) SPI (o CLK OUT 
滤波 器 ， 结 果 以 相应 的 采样 率 生 成 12 位 输出 信 TO ; l | 
号 。AD9361 主要 根据 控制 系统 发 来 的 指令 选 通 AUXADC AUXDACx XTALP XTALN 


所 需要 的 观测 频段 。 E 5 AD9361 功能 框图 
2.2 控制 系统 Fig.5 Block diagram of AD9361 

硬件 平台 采用 xilinxZynq-7000 SoC ZC702, TE 
估 套 件 ZynqSoC 采用 Xilinx 公司 28nm 工艺 的 最 新 SoC 平台 解决 方案 ZC702， 该 方案 除了 含有 Xilinx 
传统 的 高 性 能 现场 可 编程 门 阵列 外 ， 还 集成 了 两 个 ARM Cortex A9 处 理 需 内 核 ， 能 够 运行 通用 的 
Linux 操作 系统 。 本 系统 在 Linux 系统 下 能 够 方便 地 控制 快速 傅 里 叶 变 换 、 积 分 、 求 对 数 、 定 标 等 计 
算 在 现场 可 编程 门 阵列 内 部 的 运行 流程 ， 同 时 还 能 够 利用 Linux 系统 完善 的 网 络 功 能 ， 实 现 系统 与 个 
人 计算 机 、 工 作 站 等 数据 存储 设备 的 通信 。 

系统 工作 时 通过 FMC LPC 接口 向 捷 变 收发 器 AD9361 信和 号 板 发 送 控 制 信 号 ， 使 其 选 通 所 需要 的 


(D  http://www.analog.com/en/ education/education-library/ mixed, signal. dsp. design. book.html 
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观测 频段 ， 同 时 接收 AD9361 传 来 的 采样 数据 。XC7Z020 通过 其 中 的 现场 可 编程 门 阵列 部 分 实现 
AD9361 基带 信号 的 接收 ， 以 及 后 续 的 快速 傅 里 叶 变换 、 积 分 、 求 对 数 、 定 标 等 运算 ， 而 整个 数据 处 
理 过 程 在 ARM A9 运行 的 Linux 执行 程序 的 监控 下 完成 ，Linux 执行 程序 还 能 通过 适当 的 人 机 接口 ， 
实现 运算 参数 的 有 效 设 定 和 和 更改。 经 过 现场 可 编程 门 阵列 处 理 信 号 的 数据 带宽 将 极 大 地 降低 ， 达 到 
ARM A9 能 够 实时 存储 和 传输 的 水 平 ， 此 时 再 用 Linux 提供 的 完善 的 文件 系统 及 网 络 功能 实现 处 理 后 
数据 的 实时 存储 和 传输 。 


3 ”实现 方法 与 结 来 


基于 AD9361 和 Xilinx ZC702 FPGA 开发 板 的 捷 变 收发 需 系 统 ， 需 要 对 Xilinx ZC702 FPGA 开发 板 
进行 镜像 烧 写 、 配 置 和 参数 修改 ,使 其 能 够 针对 射电 源 目标 频率 进行 观测 。 
3.1 系统 构建 准备 

(D 前 端 平台 ， 将 射电 望远镜 接收 机 输出 (6 GHz 以 上 采集 中 频 信 号 ，6 GHz 以 下 可 直接 采集 ) 与 
AD9361 的 接收 端 相连 ， 试 观测 中 ， 接 于 上 海天 文 台 65m 射电 望远镜 C 波段 中 频 ，Xilinx ZC702 FPGA 
开发 板 与 AD9361 的 AD-FMCOMMS2-EBZ 板 卡通 过 FMC LPC 总 线 相连 ， 

(2) 控 制 计 算 机 ， 基 于 Windows7 64 位 操作 系统 的 个 人 计算 机 ， 用 于 上 位 机 软件 集成 及 现场 可 编 
程 门 阵列 控制 ，; 

(DSD 卡 ， 用 于 烧 写 系统 镜像 ， 容 量 8G 及 以 上 , 在 ADI 公 司 官方 网 站 的 相应 页 面 下 载 镜像 文件 
压缩 包 ， 选 择 2014_R2-2015_02_06.img.xz 下 载 后 将 其 解压 缩 ; 采用 Win32DiskImager 软件 将 解压 后 的 
镜像 文件 烧 写 入 准备 好 的 SD 卡 中 ;有 以 下 注意 事项 ， 

CD 确保 SD 卡 驱 动 需 被 选择 为 写 人 设备 ; 

D 导入 正确 的 镜像 文件 路 径 ; 

O 点 击 “write” 按 钮 写 人 镜像 文件 ; 

D 写 入 完成 后 安全 退出 SD F; 

© Z, SD 卡 中 的 镜像 还 不 能 用 于 捷 变 收发 絮 设备 的 启动 ; 

(4)SD 卡 中 包含 了 各 种 开发 板 的 启动 镜像 ， 要 在 其 中 选择 与 Xilinx ZC702 FPGA 开发 板 和 AD- 
FMCOMMS2-EBZ 板 卡 相配 套 的 镜像 文件 ， 打 开 “zynq-zc702-adv7511-ad9361-fmcomms2-3” 文 件 夹 ， 将 
其 中 的 镜像 文件 从 子 目 录 复 制 到 SD 卡 的 根 目录 下 ; 镜像 烧 写 完成 。 

3.2 启动 并 配置 IP 地 址 

(1) 接 入 相关 外 设 ， 同 时 将 AD-FMCOMMS2-EBZ 板 卡 与 现场 可 编程 门 阵列 开发 板 的 夹层 卡 口 连接 ; 

(2) 然 后 将 射电 望远镜 中 频 接 收 机 输出 与 AD-FMCOMMS2-EBZ 板 卡 RX2A 相连 ; 

(3) 在 Menu 菜单 中 打开 settings 页 面 ， 进 入 Network Connections, ， 修 改 IPv4 Settings ， 将 其 中 的 IP 
地 址 设置 为 169. 254. 188. 133 。 

3.3 应 用 程序 参数 设置 

( 1) IP 地 址 的 设置 为 现场 可 编程 门 阵列 开发 板 与 PC 通信 奠定 了 基础 ， 然 后 安装 捷 变 收发 器 应 用 
程序 IIO Oscilloscope Windows 版 至 上 位 机 中 ; 打开 应 用 程序 并 建立 连接 ,在 IP addr; port 中 输入 之 前 
配置 好 的 IP 地 址 169. 254. 188. 133 ， 选 择 确 定 进 入 AD-FMCOMMS2-EBZ 参数 配置 界面 ; 

(2) 打 开 FMComms2/3/4 选项 卡 并 运行 New Plot, 人 参数 设置 与 频 域 峰值 如 网 6， 对 比 商用 频谱 仪 得 
到 图 7 的 结果 ， 可 以 清楚 地 看 到 在 751. 4 MHz 附近 出 现 知 干 峰值 ， 捷 变 收 发 需 谱 线 观测 终端 表现 良好 。 


4 8s 论 


基于 Xilinx 现场 可 编程 门 阵列 与 AD9361 的 捷 变 收发 带 系 统 实 现 了 射电 天 文 厘 米 - 分 米 波 段 罕 带 
详 线 观测 ， 采 用 集成 化 必 片 可 以 成 功 取代 目前 传统 的 分 立 式 模拟 下 变频 天 、 模 拟 滤 波 右 、 高 速 模 数 转 
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换 融 、 高 速 现场 可 编程 门 阵列 及 高 性 能 计算 机 的 观测 系统 结构 ， 并 且 可 以 保证 天 文 观测 的 可 靠 性 、 稳 
定性 需求 。 可 以 预见 ， 这 套 系 统 将 在 未 来 的 天 文 观 测 领 域 发 挥 更 大 的 作用 。 


voltage3 


ANALOG ` 
DEVICES -= 


图 6 系统 运行 与 结果 显示 界面 
Fig. 6 System running and result display interface 


RBW 10 kHz 
Ref -53.00 dBm SWT 34s 770.007900000 MHz 


CF 770.0079 MHz 
Date: 23. JUL. 2015 13:57:31 


图 7 商用 频谱 仪 FSL26 分 析 得 到 的 数据 
Fig.7 Data analyzed by a commercial spectrometer FSL26 
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Abstract; Spectral line observation plays an important role in astronomy. High spectral resolution can 


bring much convenience to observers in exploring planet details. But for high-frequency observation, it is very 


difficult for the traditional sampling method to achieve high spectral resolution. In this article, a new method 


has been proposed. Using high performance integrated circuit to observe radio astronomy signal in centimeter- 


decimeter spectral line, this method overcomes those disadvantages of its traditional counterpart and achieves a 


higher spectral resolution. with a lower sampling bandwidth. It consists of a complete set of collection, 


processing and storage solutions. Firstly, high-frequency signal is processed by the analog down converter, 


then the low-frequency signal is transmitted to the upper computer software for FFT and integral processing; 


finally, the signal 1s saved in file format as required. In this platform, parameters for data collection are given 


by the observer through upper computer software. 


Key words: FPGA; Agile transceiver; Radio astronomy; Spectral line observation; Spectral resolution 


